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Abstract 

We embody the density matrix renormalization group method for the 19- 
vertex model on a square lattice; the 19-vertex model is regarded to be equiv- 
alent to the XY model for small interaction. The transfer matrix of the 
19-vertex model is classified by the total number of arrows incoming into one 
layer of the lattice. By using this property, we reduce the dimension of the 
transfer matrix appearing in the density matrix renormalizaion group method 
and obtain a very nice value of the conformal anomaly which are consistent 
with the value at the Berezinskii-Kosterlitz-Thouless transition point. 
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The density matrix renormalization group method (DMRG) has been used for evaluation 
of eigenvalues of Hamiltonian matrix for a one dimensional quantum system This method 
enables us to increase the system size, within the matrix size which can be handled by recent 
computer resources such as the memory size and the cpu speed. Recently the DMRG method 
has been applied to the transfer matrix of classical spin system 0. 

In the both cases of quantum Hamiltonian and classical transfer matrix, spin variables 
have discrete degree of freedom. Therefore one can make matrices with finite dimension in 
the both cases. On the other hand, a spin system with continuous degree of freedom such 
as a classical XY model provides a transfer matrix whose dimension is infinite. Fortunately 
it is known that the XY model on a square lattice A is translated to a 19-vertex model for 
which the transfer matrix can be described by discrete variables expressed by arrow variables 
@. The purpose of the present study is to embody the DMRG method for the 19-vertex 
model. We show how to reduce the dimension of matrices, which are diagonalized in the 
DMRG method, by using the ice rule of the 19-vertex model 0]. We obtain very nice value 
of conformal anomaly, c = 1.006(1), at the Berezinskii-Kosterlitz-Thouless (BKT) transition 
point by use of our method. 

The 19 kinds of arrow configurations i.e. vertices are permitted by an ice rule which is 
generalized to include no arrow on a bond. A vertex weight W(i)i) depends on the kind of 
the vertex Vi G {1, 2, ■ • • , 19} at site i. The value of Vi is determined by configuration of four 
arrows as follows: 

v i = Vi(a i} Pi,'y i ,5i), (1) 

where oij,/3j,7j and 5i denote arrows on bonds surrounding the site %. Let us express an 
up and a right arrow by +1, a down and a left arrow by —1, no arrow by 0. For instance, 
we have Vi(— 1, 0, 0, +1) = 1 as shown in Fig.l 

of ref. [Q]. The ice rule is expressed by using this expression as follows: 

at -f3i- 7< + Si = 0. (2) 
Using the weights W{v{) for 19 vertices, we can describe the partition function Z as follows: 
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^=E / n^)> (3) 

{vi} ieA 

where the summation is taken over all permitted configurations of the vertices on the lattice 
A. 

Because there exists the ice rule for the 19-vertex model, the whole transfer matrix 
becomes a block diagonal which is classified by the number of arrows incoming to one layer 
of the lattice. By using this property, we can reduce the amount of calculation in the DMRG 
method. 

We now apply the infinite DMRG method to the 19-vertex model Q . In addition to the 
vertex weight W(vi), we use a renormalized weig ht W( r \vl r) ), where means a renormal- 
ized vertex and r the number of a renormalization. As an initial value, W^(vl°^) is equal 
to W(vi). The transfer matrix for N, the number of arrows, is composed as 

tP(t)i, (3 2 , Tfe, At|6, £3, 

=J2 W^Xv^W^W^iv^Wivt), (4) 

where the number N is obtained by 

N = N r ( Vl )+/3 2 + N r ( m )+(3 4 

= N r (^) + 5 2 + N r (^) + 6 4 . (5) 

Here the number of arrows included in the renormalized vertex £j is denoted by N r (£i) and 
is equal to 8^ in the initial step of the DMRG procedure. We denote an eigenvector of 
this transfer matrix by iPn\(t)i, (3 2 , Pa) which corresponds to the k-th eigenvalue. As the 

(r) 

infinite DMRG method, we construct the density matrix pN r ( m )+i3 2 fc as f°U° ws: 

PN ) r ( m )+p 2 ,k(Vup2\£l,5 2 ) 

=>] WN,k(VU fa, m, Pd^N,k(^ S 2i Va, (6) 

Notice that the density matrix is labeled by N r (r)i)+/3 2 not by N. From Eq.(|^), since rj 3 = £ 3 
and fa = 64, N r (r]i) + (3 2 = N r (£i) + 5 2 , therefore the density matrix has a block diagonal 
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structure classified by the total number of arrows for the half system. In order to construct 
the renormalized vertex weight W^ r+l \ we diagonalize the density matrix PN r (rn)+foA anc ^ 
obtain its eigenvectors V^^+farf ■ m the present study, we use the eigenvector of the 
transfer matrix for the largest eigenvalue k — 1. The renormalized vertex state is labeled 
by rj[ which means that Vjv^j+ft,^ corresponds to the 7^-th eigenvalue of PN r (m)+/32,i- 
We determine the upper limit I of i][ as follows: 

[ gr+2 /gr+2 < \ 

^ (7) 
m (3 r+2 > m), 

where m is the number of states of the density matrix which we take into account. The last 
step of the DMRG for the 19-vertex model is a construction of the renormalized weight for 
the vertex as follows: 



^ +1) ri r+1) («i^i,a 3 ,eD) 



= J2 VN r ( v [),r,[(vi,P2)W( r \v[ r) (ai,rix,a2,(,i)) 

xW(v 2 (a 2 , fa, a 3 , 5 2 ))V Nr (ft)ft(Zi, S 2 ). (8) 

The total number of arrows included in the renormarized vertex becomes 

N r (7il)=N r (Q=N r (Z 1 ) + 8 2 . (9) 

We need W {r+1 \vi r+1) (a 3 , rf 3 , a 4 , Q) to return to the first step of the DMRG method, but 
we do not need to calculate it. Since in our method for the 19-vertex model the system has 
a symmetry of translation, we can use W^ r+1 \v[ 7+1 ^) as W^ r+1 ^(v^ +1 ^). Then we return to 
the first step represented by Eq.(£|) to iterate the DMRG procedure. By iterating the above 
DMRG procedure for the 19-vertex model with the restriction of the total number of arrows 
as mentioned here, we can make the system size L increase systematically. The advantage 
of the present method is that the dimension of the transfer matrix decreases dramatically 
by considering the conservation law of the number of arrows in this DMRG method. For 
example in the case of iV = 0,m = 35 and L — 12, the dimension of matrix is dramatically 
reduced from 531441 down to 1545. 



The value of the conformal anomaly has to be 1 for the BKT transition point. The result 
at K = 1.0866 is shown in Fig.l. The value of interaction is equal to that estimated as the 
critical value in ref. ||. We obtain c = 1.006(1) which is consistent with the value for the 
BKT point. 

As summary, we have succeeded to reduce the matrix dimension to be diagonalized in 
the DMRG method for the 19-vertex model by using the ice rule. A very nice value of the 
conformal anomaly is obtained at BKT transition point by the present approach. 
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Figure caption 
Fig.l The size dependence of the free energy. 
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